Inertial effects in Biittiker-Landauer Motor and Refrigerator at the Overdamped Limit 
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We investigate the energetics of a Brownian motor driven by position dependent temperature, 
commonly known as the Biittiker-Landauer motor. Overdamped models (M = 0) predict that 
the motor can attain Carnot efficiency. However, the overdamped limit (M —* 0), contradicts the 
previous prediction due to the kinetic energy contribution to the heat transfer. Using molecular 
dynamics simulation and numerical solution of the inertial Langevin equation, we confirm that the 
motor can never achieve Carnot efficiency and verify that the heat flow via kinetic energy diverges as 
M _1//2 in the overdamped limit. The reciprocal process of the motor, namely the Biittiker-Landauer 
refrigerator is also examined. In this case, the overdamped approach succeeds in predicting the 
heat transfer only when there is no temperature gradient. Its found that the Onsager symmetry 
between the motor and refrigerator does not suffer from the singular behavior of the kinetic energy 
contribution. 

PACS numbers: 05.40.-a, 05.40.Jc, 05.10.Gg, 05.70.Ln 



o 



in 
c3 



i 

C 

O 
o 



> 

O 

(N 

o 

oo 
o 



X 



I. INTRODUCTION 



Since the industrial revolution, thermodynamics has 
been a guiding principle for the development of new tech- 
nology. We are now entering the era of nanotechnology 
where we can construct and manipulate nanoscale ob- 
jects as we desire. Realization of nanoscale machines is 
within our reach. However, we still have to overcome var- 
ious issues. As the size of a system approaches to that of 
molecules, thermal fluctuations begin to play a significant 
role. It would be rather difficult to operate a nanoma- 
chine against strong thermal fluctuations. Instead, the 
nanomachine must be able to work harmoniously or even 
collaboratively with the fluctuations. In order to design 
such machines, we need to understand thermodynamics 
of small systems taking into account large thermal fluctu- 
ations. Furthermore, the molecular machinery in biolog- 
ical systems such as motor proteins are similarly subject 
to large thermal fluctuations. Thermodynamics, at the 
macromolecular level, is also essential in the investigation 
of such biological machines. 

Unfortunately, since it was originally developed for 
macroscopic systems where fluctuations are negligi- 
ble, standard thermodynamics is often powerless in 
cases where fluctuations dominate. We often resort to 
stochastic approaches such as Fokker-Planck equation 
or Langevin equation. Despite the fact that these ap- 
proaches have been successfully used for many years, the 
relation between thermodynamics and stochastic meth- 
ods is not well established. It was only 10 years ago 
when a general theory of energetics such as heat within 
the stochastic regime (stochastic energ etks) mug 

was 

developed. The validity of the theory needs to be system- 
atically tested with experiment or first principles simula- 
tion. 



'Electronic address: kawai@uab.edu 



When traditional thermodynamics was developed, 
Carnot engine played a key role as an idealized model. 
Similarly, Brownian motors [ij have been basic work- 
ing models for systems dominated by thermal fluctua- 
tions. In particular, autonomous thermal engines such as 
the Feynman-Smoluchowski (FS) motor [f| and Biittiker- 
Landauer (BL) motor 0, 0] , unlike other Brownian mo- 
tors, do not require time-dependent external influence. 
These motors are driven solely by thermal fluctuations 
and their motility disappears when the motors become 
macroscopic in size. 

There is no difficulty in the investigation of their motil- 
ity using standard stochastic approaches. However, it 
was not straightforward to investigate the thermodynam- 
ics of these systems. In his celebrated textbook Feyn- 
man attempted to investigate the thermodynamics of the 
FS motor and concluded that it can reach Carnot effi- 
ciency. Yet, he overlooked the effect of fluctuations and 
later it was shown that the heat transfer between two 
heat reservoirs never ceases even when the motor moves 
quasistatically and thus it is not possible to reach the 
Carnot efficiency [3, H, Q. Even when the average ve- 
locity of the motor is zero, fluctuations around the mean 
value can transport some energy. It is this fluctuation 
that transports heat. Recently, a simpler model of the 
FS motor was developed and their result confirms the 
presence of such heat transfer Similar heat transfer 
was investigated in the problems of adiabatic piston 
and shared piston [T3|. The reverse process of the FS 
motor, namely the FS refrigerator has also been studied 
using various models [H, 0, . 

The BL motor is just an overdamped Brownian par- 
ticle in a periodic potential field subject to spatially in- 
homogeneous temperature. When temperature changes 
across a potential barrier, the Brownian particle jumps 
over the barrier more often from the hot side to the cold 
side than the other way [f| 0, Therefore, Brown- 

ian particles move in one direction on average. Whereas 
the FS motor is simultaneously in contact with two heat 
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baths, the BL motor moves from one heat bath to an- 
other by itself. Therefore, the BL motor is a different 
class of Brownian motor from the FS motor. In order for 
the BL motor to operate continuously, it must be ther- 
malized with the local environment before entering the 
next heat bath. Thus, it works better in the overdamped 
limit but fails in the underdamped limit 

Like the FS motor, there is no difficulty in explaining 
the motility of the BL motor. However, the thermody- 
namics of this system is not straightforward. Intuitively 
we expect non-vanishing heat transfer even when the av- 
erage current is zero since the Brownian particles can 
move back and forth between the hot and cold regions 
by thermal fluctuations [HI, [l9[ . Despite this anticipa- 
tion, some previous investigation claimed that the BL 
motor can reach the Carnot efficiency [2Cj HH and under 
certain conditions can act as a refrigerator, attaining the 
corresponding Carnot coefficient of performance [2l| . It 
turns out, however, that the overdamped Langevin ap- 
proach fails to predict such heat transfer for the BL motor 
while it worked fine for the FS motor [l|, [i| . It appeared 
that even when the system is in the overdamped regime, 
inertial mass (M) apparently plays a critical role in cer- 
tain thermodynamic processes. In fact, heat evaluated 
by assuming M — at the beginning does not agree with 
the result obtained by taking the limit M — > at the 
end. This singularity has been phenomenologically pre- 
dicted [ljlGjl but not yet experimentally tested. Due to 
the lack of experimental confirmation, this issue is still a 
subject of debate HH, [Hj]. 

In this paper, we investigate the BL motor and its re- 
ciprocal process, the BL refrigerator, using stochastic en- 
ergetics based on Langevin approach, with and without 
inertial effects and compare the results with molecular 
dynamics simulation. Our main objective is to check the 
failure of the overdamped Langevin method and the va- 
lidity of the inertial Langevin method by comparing the 
results with molecular dynamics simulation. The over- 
damped method appears to fail for the BL motor involv- 
ing inhomogeneous temperature but works fine for the 
BL refrigerator operated with homogenous temperature. 
On the other hand, these two processes are related to 
one another through Onsager's symmetry, which is not 
obvious for a system with inhomogeneous temperature, 
prompting us to investigate the validity of the Onsager 
symmetry in the overdamped regime [231 ]. We will check 
the validity numerically, taking into account inertial ef- 
fects. 

In the next section, we introduce a concrete model of 
BL motor/refrigerator which can be investigated by both 
Langevin and molecular dynamics simulations. In sec- 
tion IIII1 a brief heuristic discussion based on the over- 
damped Langevin equation is given. In section IIVI the 
results of the overdamped approach are compared with 
numerical simulations of the Langevin equation taking 
inertia into account and the molecular dynamics simu- 
lation. Further discussion and conclusion follow in the 
final section. 
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FIG. 1: (Color online) Two rectangular reservoirs filled with 
gas particles at temperatures Ti and T2 are alternately con- 
nected. Gas particles (red and blue particles) are confined in 
the cell and only Brownian particles (large black particle) can 
move through the walls. Brownian particles are subjected to 
a piece- wise linear potential as shown in the bottom. 



II. THE MODEL AND METHODS 

A. The model 

We consider a chain of two-dimensional cells aligned in 
the x direction as illustrated in Fig. [T] Each cell has a 
width of L/2 and is filled with TV gas particles confined in 
the cell they belong to. No direct heat exchange between 
cells through the walls is permitted. The gas particles in 
a cell act as a heat reservoir with its own temperature, 
independent of the temperature in other cells. Brownian 
particles of mass M are placed in the cells. Unlike the 
gas particles, the Brownian particles are allowed to move 
freely from one cell to another through the walls and they 
are also subject to a potential field U{x), which is peri- 
odic in the x direction with a periodicity L and constant 
in the y direction. For simplicity, we use a piecewise- 
lincar potential: 



U(x) 




for < x < , 

L 

2 



for < x < L, 



(1) 



where Uq is the potential height. The x coordinate is cho- 
sen such that the location of potential maxima/minima 
coincide with the cell boundaries as shown in Fig. [TJ In 
addition to the periodic potential, a constant external 
force F is exerted on the Brownian particles. We fur- 
ther assume that temperature is periodic with the same 
periodicity as the potential and piecewise constant: 



T(x) 



Ti forO<x<-|, 
T 2 for k < x < L. 



(2) 



Temperature is measured in energy unit (ks = !)■ 
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For technical simplicity, we consider only two cells (cell 
1 for < x < L/2 and cell 2 for L/2 < x < L) with 
a periodic boundary condition instead of infinitely long 
chains. 



B. Langevin approaches 

Molecular Dynamics (MD) simulation is a useful tool 
to investigate this kind of model. However, a simpler 
mathematical model is desirable since MD simulation is 
computationally quite demanding. Apart from long-time 
fluid dynamical effects in two dimensions, the motion of 
the Brownian particle in the x direction can be investi- 
gated by the one-dimensional Langevin equation: 



Mv = -j(x)v - U'(x) + F+ y / 2j(x)T(x)£(t), 



(3) 



where x and v are position and velocity of the Brownian 
particle and £(i) is a standard Gaussian white noise: 



«(<)>= 0, =S(t-a). 



(4) 



Here, and later on, an overdot refers to derivative taken 
with respect to time and a prime means derivative taken 
with respect to space. The position-dependent friction 
coefficient j(x) is assumed to be periodic and piece-wise 
constant in the same way as temperature: 



the details of each system. For our model (overdamped 
Brownian particle subject to inhomogeneous tempera- 
ture) the appropriate Fokker-Plank equation is given by, 



dP(x,t) 
di 



d_ 

dx 



l 



(U'(x) 



~ F ^Tx T ^ 



P(x,t) 



(7) 

While solving Eq. ([7]) for the piecewise constant temper- 
ature ©, the proper boundary conditions are T\P{$) = 
T 2 P(L - S) and T Y P{L/2- S) = T 2 P(L/2 + 5), where 5 
is infinitesimally small [2a ]. 

Corresponding to the Fokker-Planck equation ([7]), the 
correct form of the overdamped Langevin equation, in 
the Stratonovich interpretation, is 



j(x)x 



-U'{x)+F 
1 d 



^2T(x) 1 {x)i(t) 



27 (x) dx 



T(x)i{x)] 



(8) 



as derived in [26|, [27|. In this paper, we shall call Eq. ([8|) 
the overdamped Langevin equation and refer to Eq. (|3|) 
as the inertial Langevin equation. 

As an indicator of overdamping, we introduce a dimen- 
sionless frictional coefficient: 



7 = 



7*o 



7 L 



(9) 



When 7 >> 1, the system is in the overdamped regime 



7(2;) 



71 for < x < j, 

72 for I < x < L. 



(5) 



It has been shown that the Langevin Eq. ([3]) correctly 
predicts the behavior of Brownian particles. However, 
due to its mathematical difficulty, further approximation 
is often used. When the relaxation time r = M/7 is much 
smaller than a mechanical time scale to = -J ML 2 /Uq, 
the inertial term in Langevin equation ^ is usually ne- 
glected. Setting M = in Eq. © we obtain a popular 
overdamped Langevin equation, 



j(x)x = -U'(x) +F+ y/2-y(x)T(x)Z(t) 



(6) 



While this equation is widely used in many different sub- 
ject areas, we must be careful with M = since both r 
and t are zero at the same time. Strictly speaking, the 
overdamped condition r <C t should be satisfied only in 
the sense of the limit M — > 0. 

Further complication arises when temperature or fric- 
tion coefficient depends on the position. Simple omission 
of the inertial term does not lead to the correct over- 
damped Langevin equation. One can also study stochas- 
tic processes in the overdamped limit by the equivalent 
Fokker-Planck equation. There is, however, no univer- 
sal Fokker-Planck equation that describes a system with 
non-uniform temperature. Van Kampen f24| found that 
a particular form of Fokker-Planck equation depends on 



C. Stochastic energetics 

We investigate the thermodynamic behavior of the BL 
system using stochastic energetics introduced by Seki- 
moto 0, [3, Heat flux from the gas particles in the 
«-th cell to the Brownian particles is defined as 



(10) 



where (■■■){ indicates ensemble average taken while the 
Brownian particles are located in the i-th cell. Using the 
inertial Langevin equation , Eq. (fTU)) can be evaluated 
in three terms: 



Qi = ^j t {i 2 ) i + {U'{x)x) i -F{x) i 

= Qf E + Qf E + Q'l 



(ii) 



where the first two terms on the rhs are the kinetic en- 
ergy and potential energy contribution to the heat flux, 
respectively, and the last term is the Joule heat. 

Work done on a Brownian particle by the external force 
in each cell is given by, 



Wi = F(x) 



(12) 



In the steady state, the net energy flux to the Brownian 
particles must be zero, and thus the energy gained by 
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FIG. 2: Mean velocity as a function of the temperature differ- 
ence AT = Ti — T2 between the two cells. The average temper- 
ature of the whole system is fixed to T av — (Ti + Ii)/2 = 0.5. 
Solid and dashed lines correspond to data obtained from the 
inertial Langevin equation ((3]) and the overdamped Langevin 
equation ([8j respectively, while circles correspond to MD sim- 
ulation. The parameter values are M/m — 4.0, as = 4.0 and 
F = 0. The dimensionless frictional coefficients vary from 
71 = 72 = 17.7 at AT = to 71 = 24.4 and 72 = 5.6 at 
AT = 0.9. As 72 becomes less overdamped with the in- 
crease of AT, the overdamped model deviates from the in- 
ertial Langevin equation. 

the Brownian particles in cell 1 must be cancelled by the 
energy loss in cell 2. Therefore, heat flux from cell 1 to 
cell 2 via the Brownian particles is defined as 

Qi^2 = Qi + W t = Qf E + Qf E . (13) 



D. Molecular Dynamics simulation 

In order to check the validity of the Langevin approach, 
we performed intensive molecular dynamics simulation. 
In our MD simulation the heat bath consists of two- 
dimensional hard disks of mass m and diameter a. The 
Brownian particle is of mass M and diameter as ■ Inclu- 
sion of external forces significantly reduces the numeri- 
cal advantage of hard disk MD simulation. However our 
algorithm (see Appendix [X| is fast enough to realize a 
sufficient number of trajectories for ensemble average. 

In order to compare the results of Langevin approaches 
with that of MD simulation, we need to find the corre- 
sponding frictional coefficient. We use an analytical ex- 



pression obtained for an ideal gas [28| : 



7 = CTBpv27rmT , (14) 

where p is the density of the gas. Note that this ideal 
7 does not depend on M. Alternately, we could use 7 
measured in MD simulation. Infact, we can get a bet- 
ter agreement between MD simulation and Langevin ap- 
proaches if the measured values are used. However, the 
measured values depend on M, hindering the real mass 
dependency of the Langevin equation. Thcrfore, we use 
the theoretical frictional coefficient (fl4]) in the Langevin 
equations (|3J) and (J5J. 

In all our MD simulations we use N = 1000 gas parti- 
cles (m = 1, a = 1) in each cell, with the size of each cell 
being 250 x 400 so that the density p = 0.01. The spa- 
tial period of the piece-wise linear potential is L — 500 
and the barrier height Uq = 1.0. As the temperature of 
the reservoirs change during the course of the simulation 
there is no well defined stationary state. However, all 
our simulations were carried out using a sufficiently large 
number of gas particles so that the system remains in a 
quasi-steady state for a sufficiently long period of time, 
enabling us to measure the various physical quantities 
properly. 



III. HEURISTIC DISCUSSION WITH THE 
OVERDAMPED MODEL 

In this section we first review the known properties of 
the BL motor and its reciprocal process BL refrigerator, 
in the overdamped limit. When the temperature of the 
cells are different (Ti > T2), the Brownian particles in the 
high temperature cell can reach a higher potential energy 
region than those in the low temperature cell. Hence, the 
Brownian particles tend to move from the hot cell to the 
cold cell over the potential barrier. At the other cell 
boundary the potential is minimum and both cold and 
hot Brownian particles can easily cross to the other side. 
Therfore, the Brownian particles flow in the positive x 
direction on average even in the absence of external force 
(F = 0). When an external load is applied (F < 0), the 
Brownian particles can do work against it as a motor. 

Asymmetry in the spatial distribution of the Brownian 
particles due to the temperature difference is the main 
driving of this motor. The overdamped Langevin equa- 
tion ijHJ) or Fokker-Planck equation ([7]) provide an ana- 
lytical expression of the spatial distribution, from which 
we obtain the average velocity of the motor (See Ap- 
pendix [B]) 
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2sinh(0i + fa>) 



sinh (fa + fa ) 



(15) 



r 




FIG. 3: (Color online) Upper panel: Steady state spatial dis- 
tribution P(x) of Brownian particles as obtained from Fokker- 
Planck equation [black line], overdamped Langevin equa- 
tion (Jsj) [red line], inertial Langevin equation ((3} [green line] 
and MD simulation [blue line]. The parameter values are 
M/m = 5.0, <T B = 6.0, 2i = 0.6, T 2 = 0.4 and F = 0. 
71 = 26.0 and 72 = 21.3 correspond to the overdamped 
regime. The two insets show the details of P(x) near the 
temperature boundaries at x/L — 0.5 and x/L — 1.0. Lower 
panel: The locally averaged kinetic energy of the Brownian 
particle along the x axis, obtained from the inertial Langevin 
equation @- 



where fi <2 — —F ± 2U /L and fa = -/ l L/(4T l ). In the 
absence of external force (F — 0), Equation (jT5J) is always 
positive for T\ > T2 and hence the Brownian particles 
move in the positive direction as expected. In fact, par- 
ticle distribution and velocity in the overdamped regime 
agree reasonably well with the results obtained from nu- 
merical simulation of the inertial Langevin equation and 
molecular dynamics simulation (see Figs. [21 and [3]), sug- 
gesting the validity of the overdamped Langevin equa- 
tion. 



Now, we consider the thermodynamic efficiency of the 
BL motor. Based on the successful prediction of par- 
ticle distribution and velocity it is natural to use the 
overdamped Langevin approach for other quantities such 
as efficiency. Applying the overdamped Langevin equa- 
tion ijHJ) to the definition of heat (fTT)|) , the heat flux from 
the heat bath to the Brownian particles in cell 1 is given 
by 



1 



27 (x) dx 



[ 7 (x)T(x)]x> +{U'(x)x) 1 -F{x) 1 



(16) 

Comparing this equation with Eq. (jlip , the first term in 
the rhs of Eq. (|16|) corresponds to the kinetic energy con- 
tribution Q^ E . However, due to periodic boundary con- 
dition and Eq. (|14[) . this kinetic energy contribution van- 
ishes and we obtain Q\ = (2Uq/L — F)(x)\. Hence, 
we find the efficiency of the motor 



-W 



-IF 



2U /L - F ' 



(17) 



where W — F(v). When the motor is in a stalled state, 
(fa + fa = 0), we can show that the efficiency reaches 
the Carnot efficiency 77c = 1 — T2/J1. 

This result is puzzling since at the moment a Brow- 
nian particle enters the cold bath it is carrying Ti/2 of 
kinetic energy. When it is thermalized with the cold bath, 
Qi^2 (^1 ~F-i)l2 of heat dissipates into the cold bath. 
This heat dissipation takes place whenever the Brown- 
ian particle crosses a temperature boundary and it is 
irreversible. Due to diffusion the crossing occurs even 
when the average velocity vanishes, and this irreversible 
heat due to thermal fluctuation of Brownian particles 
never ceases as long as there is a steep temperature gra- 
dient. Derenyi-Astumian [l8| and Hondou-Sekimoto jl9l] 
have pointed out that the kinetic energy contribution is 
the dominant channel of heat transfer between two cells 
and thus the efficiency of the BL motor cannot reach the 
Carnot efficiency. This situation is similar to the case of 
the FS motor 

An interesting question is why the overdamped 
Langevin approach failed for the BL motor whereas it 
worked fine for the FS motor. Derenyi-Astumian [l8| 
and also Hondou-Sekimoto [19] phenomcnologically an- 
swered this question. Consider the thermal velocity v t h = 
y/T/M. When a Brownian particle crosses the temper- 
ature boundary, it will not be immediately thermalized 
and there is a narrow region £ t fo = v t hJ = \TM /j where 
the average kinetic energy of the Brownian particle does 
not coincide with the temperature of the reservoir (see 
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lower panel of Fig. [3]). In the overdamped limit (M — > 0) 
this region vanishes, justifying the use of the overdamped 
approach. However, the effective temperature gradient is 
\T\— T 2 \jlth oc A/ -1 / 2 . Hence, the heat current is propor- 
tional to M -1 / 2 and diverges as M — > 0. This singularity 
implies that the overdamped regime that assumes M = 
is not equivalent to the overdamped limit M — > 0. Ap- 
parently, the FS motor does not have this problem 
One of our objectives is to verify this singular behavior 
by numerical simulations. 

Next, we turn to the BL refrigerator. We assume T\ = 
T2 = T and apply a weak external force < F < 2Uo/L 
to the Brownian particle. Due to the external force, the 
Brownian particle drifts in the direction of F and its av- 
erage velocity in the overdamped limit can be obtained 
from Eq. (|15p. Unlike the motor case, the temperature 
is uniform throughout the system and thus the kinetic 
energy does not contribute significantly to the heat ex- 
change between the cells. Again using the overdamped 
model, the heat flux absorbed by the Brownian particle 
from the i-th reservoir is obtained from Eq. (|16p as 



OTT OTT 



F(x) 



(18) 

where the first term is the potential energy contribu- 
tion Qf E and the second term the Joule heat Qj. Since 
(jf E oc F and Qj oc F 2 we can have positive Qi for suf- 
ficiently small F, indicating that cell 1 is refrigerated. 
Since there is no singular behavior due to temperature 
change, the overdamped model may be good enough. 
However, as soon as refrigeration induces a tempera- 
ture difference between the cells, the heat transfer due 
to the kinetic energy reduces the power of refrigeration. 
Our second objective is to investigate whether the over- 
damped model is sufficient and if the refrigeration can be 
sustained against the heat leak due to the kinetic energy. 

The BL motor and refrigerator are a result of cross 
effect between the external force F and the temperature 
difference AT. Assuming they are small, we can make 
a connection between the motor and refrigerator using 
linear irreversible thermodynamics expressed by 



/ V r F r AT 

W = L n y + £12-^2", 

F AT 
Qi^2 = L 2 i— + £22-^5-, 

where the transport coefficients are defined as 



(19) 



L n =T 



d(v) 



T12 = T- 



d(v) 

a at' 



J 21 



T 



dF 

dF ' " <9AT 



dQi^ 2 T 2 <9Qi^2 
1^22 — 1 



(20) 



Here, the partial derivatives are evaluated at AT = 
and F = 0. In the overdamped case we can obtain all 
Onsagcr coefficients from Eqs. (fTB"|). (TT51) . and (TTB]), as 



C/2 

Lu = -7jt£, L 2 2 - 



U 4 

L \2 



TL ' 



(21) 



where £ = [47sinh 2 ([/o/2T)] _1 . These coefficients sup- 
port the Onsager symmetry L12 = L21. Moreover, 
Eq. (f2Tjl indicates T11 L22 — L12L21 = 0, implying no en- 
tropy production [291 ] and thus the motor operates at the 
Carnot efficiency. Hence the overdamped model again 
erroneously predicts the highest possible efficiency. The 
singular behavior of Qi^2 at M = casts some doubt 
on the Onsager symmetry since the two limits M — > 
and AT — > do not commute. Our third objective is to 
numerically investigate the validity of linear irreversible 
theory and the Onsager symmetry at the overdamped 
limit. 



IV. RESULTS 
A. Motor 

First we show in Fig. [3] that the numerical solution of 
overdamped Langevin equation ([5]) agrees perfectly with 
the solution of the corresponding Fokkcr-Planck equa- 
tion ([7|) , confirming that the extra term in Eq. ([8]) is nec- 
essary. Furthermore, the results of overdamped Langevin 
equation are in good agreement with the numerical re- 
sults from inertial Langevin equation (J3]) and molecular 
dynamics simulation, except for narrow regions at the 
temperature boundaries (see the insets in Fig. [3]). Al- 
though it is in general small, this error in the overdamped 
regime should not be overlooked. 

As we discussed in section IIII1 the overdamped 
Langevin method assumes that the Brownian particles 
are locally in equilibrium with the the thermal reservoir. 
Therefore, the average kinetic energy of the Brownian 
particles is assumed to be the same as the local tem- 
perature of the heat bath. This assumption, however, 
fails near the temperature boundaries since the Brow- 
nian particles entering from one temperature region to 
another are not immediately thermalized with the new 
environment even when the mass is very small as Fig. [3] 
indicates. The size of the transition region is the ther- 
malization length £ t h- In the overdamped limit [M — > 0), 
indeed the transition region vanishes but only slowly as 

Vm. 

It turns out that such a subtle error in the overdamped 
Langevin method does not lead to large errors in cer- 
tain quantities obtained from it. Figure [5] shows that 
the overdamped Langevin equation predicts the average 
velocity as well as the inertial Langevin equation does. 
Due to this success we feel safe to use the overdamped 
approach to compute other quantities. However, the cal- 
culation of the heat is a different story. As we discussed in 
the previous section, the overdamped approach predicts 
that the heat transfer between the cells is only from po- 
tential energy contribution. However, Hondou-Sekimoto 
(l9| pointed out that the narrow transition region plays 
a dominant role in heat transfer and the overdamped 
Langevin equation fails in the investigation of heat cal- 
culation. The inertial Langevin approach concludes that 
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FIG. 4: (Color online) Mean velocity (a), heat flux due to po- 
tential energy (b), and heat flux due to kinetic energy (c) as 
a function of the mass ratio M/m. Circles, dashed lines and 
solid lines correspond to MD simulation, overdamped model 
and the inertial Langevin equation, respectively. The param- 
eter values are Ti = 0.6 and T 2 = 0.4, a B = 5.0 and F = 0.0. 
As the mass ratio changes from M/m — 0.1 to M/m = 60, 71 
varies from 48.5 to 6.3, and 72 varies from 39.6 to 5.1. In (c) 
the dotted line indicates 0.000215(M/m)~ 1,/2 , an empirical fit 
to the result of inertial Langevin equation. 



the kinetic energy contribution is actually dominant. 

Using scaling arguments, Derenyi-Astumian [l8 j and 
Hondou-Sekimoto [l9( predicted that the kinetic energy 
contribution diverges as M -1 / 2 at M — 0. Figure H] 
clearly shows such a singularity in good agreement with 
their theory and confirms Q^ E ^> Q^ E . We also find from 
Fig. [4] that result obtained from the inertial Langevin 
equation approaches the MD simulation result as the 
mass ratio M/m increases. Even for small M/m, the 
inertial Langevin equation predicts the correct order of 
magnitude of the motor velocity and the heat flows, 
though, in general, the Langevin approach is not ap- 
plicable in such a case. The motor velocity and heat 
flow via potential energy predicted by the overdamped 
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FIG. 5: The efficiency of the motor vs the magnitude of the 
external load, for mass ratios M/m = 5.0, 20.0 and 50.0. Solid 
lines correspond to numerical solution of the inertial Langevin 
equation and circles represent MD simulation (for M/m — 5.0 
only). Other parameter values are T\ — 0.7, T2 = 0.3, and 
<jb = 5.0. The dimensionless friction coefficients are 71 = 
23.4, 72 = 15.3 for M/m = 5.0; 71 = 11.7, 72 = 7.7 for 
M/m = 20.0; and 71 = 7.4, 72 = 4.9 for M/m = 50.0. 



model, deviate from the inertial Langevin equation and 
MD simulation result as the mass ratio increases be- 
cause the inertial effect becomes too large to be ignored. 
The good agreement between the inertial Langevin equa- 
tion and the molecular dynamics simulation, shows that 
Sckimoto's definition of stochastic energetics predicts the 
heat correctly even when temperature is spatially inho- 
mogeneous. 

In Fig. [5J we plot the efficiency of the motor normal- 
ized by Carnot efficiency 77c j asa function of the exter- 
nal load F normalized by the magnitude of the force ex- 
erted by the periodic potential energy. In contradiction 
to the overdamped model, the efficiency is far below the 
Carnot limit. While this result was phenomenologically 
argued for in [lH Il9| , here we confirmed with molecular 
dynamics simulation and numerical solution of inertial 
Langevin equation that the kinetic energy contribution 
greatly reduces the efficiency. Even when the motor is 
operated at the quasistatic limit with a stall force, the 
irreversible heat transfer via kinetic energy persists and 
thus the Carnot limit is unattainable. Furthermore, the 
divergence of Qf E at M = diminishes the efficiency of 
the motor to zero, in the overdamped limit. 



B. Refrigerator 

In the refrigerator mode, two cells initially have the 
same temperature and the Brownian particles are driven 
by an external force F. In Fig. [6l MD simulation illus- 
trates the cooling of cell 1 at the expense of the heating of 
cell 2. Note however, that the average temperature T av 
increases with time as Joule heat is dissipated in both the 
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FIG. 6: (Color online) Temperature of cells 1 and 2 and the 
average temperature T av = (Ti + T<i)l1 as a function of time 
obtained from MD simulation. Initial temperatures of cells 
1 and 2 are T\ = T2 = T = 0.5 and an external force 
FL/(2Uo) = 0.5 is applied to the Brownian particle in the 
positive direction. The other parameter values are M/m — 20 
and ob = 8.0, leading to a dimensionless frictional coefficient 
71 = 72 = 15.9 (overdamped regime). Temperature drops at 
the beginning because the Brownian particles are initially at 
rest. 
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FIG. 7: (Color online) The various components of Q\ (from 
top to bottom the heat flows are Q™, Qi, Q^ E and Q\ re- 
spectively) of the refrigerator as a function of the external 
force. The parameter values are M/m = 20.0, ob = 8.0 
and Ti — T2 = 0.5 so that 71 = 72 = 15.9 (overdamped). 
Circles represent MD data, dashed lines correspond to the 
overdamped Langevin equation and solid lines represent the 
inertial Langevin equation. 



cells. In the Langevin approach, we cannot see such tem- 
perature changes since the temperature needs to be kept 
constant. However, it allows us to investigate heat trans- 
fer between the cells. Figure [7] shows the components 
of heat from cell 1 to the Brownian particles as a func- 
tion of F. The potential energy contribution increases 
linearly with F whereas Joule heat decreases (Q'[ < 0) 
as F 2 [see Eq. (fT5|) ]. For sufficiently small F, the po- 
tential energy contribution wins and hence the Brownian 
particles extract heat from cell 1 and dump it in cell 
2. There is an optimal F at which the cooling effect 
is maximum. The agreement between the overdamped 
and inertial Langevin equations and with the molecular 
dynamics simulation is good for Q\ E and Q\. 

Since temperature is uniform, we don't expect signifi- 
cant kinetic energy contribution. However in Fig. [Jj we 
found that there is a small amount of kinetic energy con- 
tribution which is absent in the overdamped approach. 
In the overdamped regime we assume that the change 
in potential energy immediately dissipates into the reser- 
voir. However, with a finite mass, the potential energy 
is first converted to kinetic energy which dissipates at a 
later time. For example, when a Brownian particle slides 
down the potential slope near a cell boundary, the po- 
tential energy change is transported to the next cell as 
kinetic energy where it is dissipated as heat. The amount 
of "unthcrmalizcd" kinetic energy the Brownian particles 
acquire from the potential energy is approximately, 

Qf E « -U'(x)£ t M = -^^VmT. (22) 
which is linear with F through (i) in agreement with the 



result of the inertial Langevin equation (see Fig. [7]). Both 
this kinetic energy contribution and the potential energy 
contribution are in proportion to F, yet the magnitude 
of Qi is negligibly smaller than that of Q\ E at the over- 
damped limit. Therefore, the refrigeration is still possi- 
ble. Figure [8] shows that the PE contribution approaches 
the overdamped model as M — > 0. The KE contribu- 
tion approaches to zero as (M/m) 1 / 2 , in good agreement 
with the phenomenological prediction (f2"2"| . Figures [7| 
and [8] both show that the overdamped model is in good 
agreement with the inertial Langevin equation as well as 
MD simulation in predicting the velocity and heat flows. 

The overdamped Langevin equation seems to work well 
for the refrigerator. However, the cooling effect creates 
a temperature difference between the two reservoirs and 
the overdamped approach again fails. In turn, the tem- 
perature difference induces a thermodynamic force op- 
posing the external force F, through the Brownian mo- 
tor mechanism. As the temperature difference increases, 
the motor and refrigerator effects eventually cancel each 
other, so that the average velocity of the Brownian par- 
ticle becomes zero and the cooling ceases. 



C. Onsager symmetry 

As discussed in section IIII1 we need to evaluate the 
transport coefficients (f2"0"|) with a certain care. Since 
the heat transfer diverges when the limit M — > is 
taken before AT — > 0, we first evaluate the coefficients 
at the finite mass numerically from the response curves 
like Figs. [2] and [7J Then, we reduce the mass toward 
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FIG. 8: (Color online) The various components of Qi of the 
refrigerator as a function of the mass. The parameter val- 
ues are a B = 5.0, Ti = T 2 = 0.5 and FL/(2U ) = 0.375. 
The dimensionless friction coefficients (71 =72) vary from 
31.3 at M/m = 2.0 (overdamped) to 5.7 at M/m = 60.0 
(weakly damped) . Refer to Fig. [7] for an explantion of the 
various data lines and symbols. The dotted line represents a 
fit (oc y/ M/m) to the kinetic energy contribution according 
to Eq. (122)) . 

M = 0. Figure [5] shows that as M decreases, the prod- 
uct of off-diagonal coefficients approaches to the result of 
overdamped Langevin equation obtained from Eq. (|21[) . 
However, the product of diagonal coefficients diverges as 
(M/m) -1 / 2 , reflecting the divergence of Q KE . As ex- 
pected, the off-diagonal coefficients do not suffer from 
the singular inertial effect. Furthermore, Fig. [5] shows 
that L11L22 3> L12L21 for all values of M/m, indicating 
large entropy production consistent with the very low ef- 
ficiency. Finally, Fig. [10] numerically verifies the Onsager 
symmetry for all masses within the accuracy of simula- 
tion and that the overdamped model (f2Tj) predicts the 
correct overdamped limit. 



V. DISCUSSION 

In this paper, we investigated the thermodynamic 
properties of the BL motor and refrigerator using 
Langevin equations and molecular dynamics simulation. 
For mechanical properties such as the velocity of the 
Brownian particles, we observed reasonable agreement 
between overdamped Langevin equation and molecu- 
lar dynamics simulation. However, the overdamped 
Langevin equation failed to predict thermodynamic prop- 
erties such as the heat transfer. On the other hand, 
we found good agreement between the inertial Langevin 
equation and molecular dynamics simulation even in the 
overdamped regime. Therefore, we conclude that the 
inertial mass plays a significant role even in the over- 
damped limit. The main effect of the inertial mass is the 
kinetic energy contribution to the heat transfer. We con- 
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FIG. 9: (Color online) The product of the Onsager coeffi- 
cients T11T22 (red curve) and T12T21 (blue curve) as a func- 
tion of M/m obtained from the inertial Langevin equation. 
Data has been normalized with respect to the overdamped 
value. The dashed line represents Onsager coefficients ob- 
tained from overdamped results. The parameters for obtain- 
ing I/12 and L22 were AT = 0.025 and F = 0, while Lu and 
L21 were obtained at AT = and FL/(2U ) = 0.025, where 
Tx = (1/2)(1 + AT), T 2 = (1/2)(1 - AT) and a B = 5.0. We 
use a logarithmic scale for the y-axis. The dotted line is the 
phenomenological fit 18.4(M/m) _1,/2 , to the product L\\Lrz- 



firmed that the kinetic energy contribution is dominant 
when two cells have different temperature and verified 
a previous phenomenological prediction, that the irre- 
versible kinetic energy contribution diverges as Af -1 / 2 
at the overdamped limit. 

Recently, Van den Broeck [30| investigated the effi- 
ciency of Brownian motors using linear irreversible ther- 
modynamics and concluded that, in principle, Carnot ef- 
ficiency can be attained. He argues that when the load F 
is small, the temperature difference AT drives the motor 
in the forward direction and at the same time the mo- 
tor transfers heat from the high temperature reservoir 
to the low temperature reservoir. As the load exceeds 
the stall force, the motor moves backward, transferring 
heat against the temperature gradient as a heat pump. 
In an ideal system, the direction of velocity and heat are 
reversed at the same magnitude of the external force. 
Hence, the thermodynamic fluxes (v) and Qi_>2 simulta- 
neously vanish despite F and AT not being zero. This 
can only happen when LnL 2 2 — Li 2 L 2 i = and thus 
Carnot efficiency is achieved. 

In our case, the direction of Q FE coincides with the 
direction of (v) , and both Q PE and (v) vanish at the stall 
force. However, Q KE flows from the hot to cold reservoir 
regardless of the direction of (v). Therefore, Q KE not 
necessarily vanishes. In fact, it never vanishes and the 
Carnot efficiency cannot be achieved in the BL system. 
The efficiency of the BL motor is significantly reduced 
by the kinetic energy contribution and vanishes at the 
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FIG. 10: (Color online) The Onsager reciprocity coefficient 
I/21 vs L\2 for various M/m. The parameter values are the 
same as in Fig. [9] The solid circles represent data obtained 
from the inertial Langevin equation. The value obtained from 
the overdamped model is shown with a solid square. 
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APPENDIX A: HARD DISK MOLECULAR 
DYNAMICS SIMULATION 

Our hard disk molecular dynamics simulation is based 
on a usual event-driven algorithm [33j . In the case of 
collision between two gas particles the collision time can 
be obtained analytically as given in [33[. However, since 
the Brownian particle is subject to forces due to the piece- 
wise linear potential U(x) and the external load F, the 
calculation of collision time between a Brownian particle 
and a gas particle is not straightforward. Suppose that 
we know the relative position r and relative velocity v 
of a Brownian particle with respect to a gas particle at 
time tg. We want to find the time t = to + e, at which 
they collide. The distance between two particles at the 
moment of impact is d = (a + ctb)/2 and hence 



overdamped limit. 

Derenyi-Astumian pj§] have argued that if we could 
prevent the re-crossing of the Brownian particle over 
the temperature boundary, such irreversible heat transfer 
could be arbitrarily reduced and thus the efficiency would 
be greately improved. They proposed to place a gate 
at each temperature boundary which prevents Brownian 
particles from crossing the boundary back and forth mul- 
tiple times during a short period of time. If such a gate is 
possible, the motor could reach Carnot efficiency. How- 
ever, it is not clear at present how to construct such a gate 
without reducing the particle velocity. Based on a naive 
consideration one might expect significant reduction of 
the kinetic energy contribution by optimizing the poten- 
tial profile and the location of the temperature boundary. 
Unfortunately, the reduction of heat transfer always re- 
sults in the reduction of the motor velocity as well and 
the gain in the efficiency is very limited [311 ]. 

Recently, Humphrey et al. [32] ] proposed a quantum 
Brownian heat engine which achieves Carnot efficiency. 
In this model, reversible particle exchange between two 
cells is achieved by filtering the energy of the particle, 
without violating second law. The filter only allows par- 
ticles having a certain energy, for which the Fermi-Dirac 
distribution in the two cells coincide, to pass through. 
This ensures that the particle flow does not alter the 
thermal distribution in cither cell so that there is no ki- 
netic energy contribution. Hence, the heat engine attains 
Carnot efficiency in the quasistatic limit. Whether a sim- 
ilar heat engine can be devised in the classical regime is 
not clear. 



ve 



= d 



(Al) 



where a is a constant acceleration of the Brownian par- 
ticle. (Note that the gas particle has no acceleration.) 
To find e, we rewrite Eq. (|A1|) as a quartic equation for 
a 



Ae 3 + Be 2 + Ce + D = 



(A2) 



where A = 4v-a/a 2 , B = 4(a • r + v 2 )/a 2 , C = 8r ■ v/a 2 , 
and E — 4(r 2 — d 2 )/a 2 . Using a standard method [3J], 
we find four solutions 



A 

T 



±(f±Vg + h), 



A 

T 2 



where, 



F 
G 

H 



A 2 



- B + x 



-A 2 - F 2 -2B 
4 

AAB -8C-A 3 



l - (f±Vg=h 

(A3) 



(A4) 
(A5) 



AF 



1\fx 2 



4D 



if F^0 



if F = 



(A6) 



and x is a real solution of an auxiliary equation, 



Bx 2 + (AC - 4D)x + ABD -C 2 - A 2 D = . (A7) 



The analytical solutions to the cubic equation (IA7[) are 
given in [351 ]. 
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The correct collision time corresponds to the smallest 
positive real root of Eq. (|A3|1 . if it exists. Despite hav- 
ing analytical solutions this algorithm occasionally fails, 
particularly when the acceleration a is significantly small 
because some of the coefficients become extremely large 
causing bit-off errors. To overcome this difficulty, we im- 
prove the accuracy of the roots by iterating the Newton- 
Raphson [35j steps a few times starting from the ana- 
lytical solution to Eq. (|A7|) . We found that the present 
algorithm is faster than directly solving the quartic equa- 
tion (|A2[) by the Newton-Raphson method. 



where the net force is defined as f\ :2 = —F±2U~o/L. We 
consider only weak external forces |F| < F c — 2Uq/L so 
that it does not destroy the potential barrier and hence 
we assume fi ^= 0. Then, we find general solutions: 



/'if.'' > C\^P(~ 



P2{x) = C 2 exp 



f2(x-L) 



Do 



(B3) 
(B4) 



APPENDIX B: ANALYTICAL EXPRESSIONS IN 
OVERDAMPED CASE 

We derive analytical expressions for steady state den- 
sity P{x) and current J(x). The Fokker-Planck equation 
for steady states is simply dJ(x)/dx — 0, where the par- 
ticle current is defined by 

7 (x)J(x) = -(U'(x) - F)P(x) - ±[T(x)P(x)]. (Bl) 

ax 

In our model, the Fokker-Planck equation for the i-th cell 
is given by 

Pl'(x) + ^P;(x) = 0, (i = l,2), (B2) 



where C, and Di are the constants of integration. The 
corresponding currents are Ji — —Difi/^/i, (i = 1,2). 
Since temperature and frictional coefficient are discontin- 
uous at the cell boundaries, the density does not have to 
be continuous. However, the current must be continuous 
(Ji = J 2 = J). From Eq. (|B1|) . we find the magnitude 
of the density discontinuity as T 1 Pi(L/2) = T 2 P 2 (L/2) 
and TiPi(O) = T 2 P 2 (L) [25]. These boundary conditions 
along with normalization of the density are sufficient to 
determine the integral constants. After obtaining an ex- 
pression for J and the integral constants, we obtain a gen- 
eral expression for the average velocity [Eq. (|15[) ] from a 
general relation (v) = JL. The particle density in each 
cell is obtained from Eqs. (|B3|) and (|B4|) . 
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